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FPU phenomenon for generic initial data 
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The well known FPU phenomenon (lack of attainment of equipartition of the mode-energies at 
low energies, for some exceptional initial data) suggests that the FPU model does not have the 
mixing property at low energies. We give numerical indications that this is actually the case. This 
we show by computing orbits for sets of initial data of full measure, sampled out from the micro- 
canonical ensemble by standard Montecarlo techniques. Mixing is tested by looking at the decay 
of the autocorrelations of the mode-energies, and it is found that the high-frequency modes have 
autocorrelations that tend instead to positive values. Indications are given that such a nonmixing 
property survives in the thermodynamic limit. It is left as an open problem whether mixing obtains 
within time-scales much longer than the presently available ones. 

PACS numbers: Valid PACS appear here 
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By the "standard" FPU phenomenon we mean the cel- 
ebrated one observed for the first time in the year 1955 
[l| . Namely, in numerical integrations of the equations of 
motion for a chain of particles coupled by weakly nonlin- 
ear springs, equilibrium is not attained within the avail- 
able computational time, and a kind of anomalous pseu- 
doequilibrium does instead show up. This is observed 
at low energies, for initial data very far from equilib- 
rium. The quantities studied were the energies Ej{t) 
of the normal modes of the linearized system, and their 
time-averages Ejit) were found to relax each to a dif- 
ferent value rather than to a common one, against the 
equipartition principle (see especially the last figure of 
the original FPU report). 

It was later found by Izrailev and Chirikov [2] that 
the phenomenon disappears, i.e., energy equipartition is 
quickly attained, if energy is large enough. A long de- 
bate then followed [1, |3j H [1, 0] concerning the questions 
(still unanswered) whether the phenomenon persists in 
the "thermodynamic limit" (i.e., when the number N of 
particles and the energy E both grow to infinity with a 
finite value of the specific energy e = E/N), and whether 
it can be interpreted in a metastability perspective [1, Q . 
Another still open problem is whether the phenomenon 
persists when the dimensions are increased, passing from 
a chain of particles to a 2- or a 3-dimensional lattice [l^ . 

In the present letter we address a further problem, 
namely whether some analog of the FPU phenomenon 
occurs if generic initial data are taken rather than some 
very special ones (see also [ill El; [I!])- More generally, 
we would like to look at the problem of the approach 
to equilibrium from the viewpoint of ergodic theory, in 
which one considers in principle all initial data, weighted 
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with an invariant measure such as the microcanonical 
one. Now, in ergodic theory it is well known that an ap- 
proach to equilibrium is guaranteed if a system is proven 
to be mixing. Let us recall this. For a function / on phase 
space, define f{t) = f o <!>*, where $* is the flow induced 
by a given Hamiltonian. Denote also by ( • ) expectation 
with respect to the microcanonical measure. Then, mix- 
ing amounts to requiring (see [l3|, theorem 9.8) that for 
all (square-integrable) functions / and g the correlation 

C(i) = (/W5(0))-(/W)(g(0)) 

tends to zero as t oo. So, equilibrium occurs if the 
correlations between all pairs of functions are proven to 
decay to zero with increasing time. 

The original FPU results, although expressed in terms 
of time-averages and observed only for a very special 
class of initial data, suggests that, at low energies, the 
1-dimensional FPU system "does not have mixing prop- 
erties up the considered time" . In the present letter 
we give strong numerical indications that this is actu- 
ally the case, even in the thermodynamic limit. This is 
obtained by computing the correlations of suitable func- 
tions, the averaging being performed over initial data 
sampled out from the microcanonical ensemble through 
a suitable Montecarlo method. According to the compu- 
tations, for low enough energies the correlations appear 
to relax to some positive values. This we call an FPU- 
like phenomenon. Such a phenomenon seems to suggest 
a positive property, namely that the system did actually 
relax to some well defined anomalous state. We leave 
for further studies the questions whether such a result 
should be interpreted in a metastability perspective, and 
whether it persists for 2 and 3-dimensional lattices. 

For what concerns the functions to be investigated, we 
started up by following FPU, and restricted our attention 
to the normal-mode energies Ejit), i.e. we studied the 
autocorrelations 

C,{t) = {E,{t)EM) - {E,{t)){EM) ■ 
It will be shown later, however, that a major role is played 
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FIG. 1: The normalized autocorrelation function Cj(t)/Cj{0) of the mode-energies Ej versus time, for some selected values 
of j (j ~ 16fc, k = 0, . . . ,32) and N = 511 at two values of the specific energy e. Left panel, t = 3.16 10~^, right panel 
e = 3.16 10"^ 



by other related quantities, i.e., the energies Ej of "pack- 
ets" of nearby modes, whose relevance was pointed out 
in [H (see also [H]). 

The aurocorrelations Cj{t) of the mode-energies Ej 
were numerically estimated by integrating a sufficiently 
large number K of orbits (actually, K = 10000, apart 
from two cases which are mentioned later) , and comput- 
ing at any time the arithmetic mean of the values cor- 
responding to the single orbits. The single initial data 
were sampled out from a microcanonical ensemble at spe- 
cific energy e. This was actually implemented as follows. 
Each initial datum was extracted from a Gibbs ensemble 
(with the quadratic part only of the total Hamiltonian) 
at temperature e, and was then rescaled to let it fit the 
constraint H = Ne {H being now the total Hamiltonian). 

We took the standard a-FPU Hamiltonian, namely, 

N 2 ^ 

H{pi, . . .PN,xi, . . .,xn) ^y^P^ +yz v{xk+i - , 

k=l k=0 

with Xo = Xn+i = 0, where pk is the momentum conju- 
gated to the particle's position Xk, and the interparticle 
potential is V{r) = + ar'^/3. Units were so chosen 
that m = 1, and a = 1/4. The integrations were per- 
formed with the standard leap-frog (or Verlet) method, 
with tipical step 0.05. 

The analog of the FPU phenomenon (together with 
its disappearing at high energies) is exhibited in 
Fig. [U where the normalized autocorrelation functions 
Cj{t)/Cj{0) of the mode-energies Ej are plotted versus 



time, for some selected values of j, with N = 511 and 
a sample of 10000 initial data. The left and the right 
panels correspond to a case of a "high" specific energy 
and to a case of a "low" specific energy respectively, pre- 
cisely, e — 3.16 10~^ and e = 3.16 10~^. It is seen that 
in the case of a high energy all autocorrelations decay 
to zero essentially within the same characteristic time, of 
the order of 10^. In the case of a low energy, instead, 
the decay to sero occurs only for some modes (with a 
characteristic time of the same order of magnitude as in 
the previous case), whereas for the remaining modes the 
autocorrelations appear to have relaxed within that time 
to some asymptotic nonvanishing values c*{j). 

The natural question then arises of understanding 
whether there is any regularity in the distribution of the 
asymptotic values among the modes. We found the in- 
teresting result that the relevant parameter is the mode- 
number j /N (which is a monotonic increasing function 
of the corresponding frequency tOj ) . This is illustrated in 
Fig. O where the asymptotic values c* of the normalized 
autocorrelations are plotted versus j /N. The figure refers 
to the same values of e, N and K (number of initial data) 
as in the right panel of Fig. [TJ The first interesting fea- 
ture is that the data appear to lie on some smooth curve. 
Moreover, the shape of the curve shows that the low-j 
modes (i.e., the low-frequency ones) are the ones that 
exhibit a quick relaxation to the "final" expected value 
0, while the high-frequency modes remain "frozen" near 
the initial value 1. One should notice that it is precisely 
by looking at the correlations that the frequency can be 
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FIG. 2: The "asymptotic values" c*(j) of the normaUzed au- 
tocorrelation functions Cj (t) /Cj (0) versus j/N, for the same 
parameters of Fig. [U namely, A'^ = 511 and e = 3.16 10~^. 
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FIG. 3: The "asymptotic values" c*(j) of the normalized au- 
tocorrelation functions Cj{t) /Cj{0), versus j^/^/e, for TV — 
511 and e = 3.16 10"^ 1.0 10~^ 3.16 W*. 



found to play any role, because the microcanonical ex- 
pectations of the energies are instead all equal (equipar- 
tition). On the other hand, as the correlations are well 
known to play a major role in thermodynamics according 
to the fluctuation-dissipation theorem, one may conjec- 
ture that the anomalous behaviour discussed here might 
be of physical interest, for example for some phenomena 
of anomalous decay observed in recent experiments (see 

M)- 

We come now to the dependence of the function c*{j) 
on the specific energy e. The curve is expected to reduce 
to the straight lines c* = and c* = 1 for large and 
small values of e respectively. We found the interesting 
result that, for a fixed N, the curve is a function of only 
one variable, precisely, one has c*{j,e) = /(j/-\/e). This 
is illustrated in Fig. [3l where, still for N = 511, c* is 
plotted versus {j/N)/y/e for three values of e, namely, 
e = 3.16 10-3, 1.0 10-3, 3.16 lO"'^. A rather good super- 
position of the curves seems to be observed. In particular 
notice that, with increasing e, the domain of the curve 
shrinks to the left, so that c* is found to approach the 
value zero. Thus for large e one has a complete decay to 
zero of the correlations, i.e., an analog of the Izrailev- 
Chirikov phenomenon. Notice that for each e the asymp- 
totic values c*(j) had to be evaluated at a suitable time, 
namely that at which stabilization just started occurring, 
in the sense of FigU] Such a relaxation time was found 
to increase as 1/e with decreasing e. 

The last point we address concerns the dependence 
of the results on the number N of particles. This is a 



quite delicate matter, on which we feel we got an inter- 
esting result. To begin with we point out that, if one 
takes a naive approach and plots the curves analogous 
to that of Figl2] for increasing values of N, the curves 
are found to collapse towards the trivial one c* — 0. 
This is shown in Fig|31 left panel, where the curves for 
TV = 511, 1023, 2047, 4095 are reported, for the same e as 
in Fig. [21 Concerning the number K of initial data, this 
had forcedly to be diminished with increasing TV, and we 
had to pass from K = 10000 for TV = 511 and 1023 to 
K = 5000 and 2000 for TV = 2047 and 4095 respectively. 
This, by the way, explains the broadening of the "curves" 
for the two large values of TV. 

It would however be incorrect to infer from such a col- 
lapse that mixing occurs in the thermodynamic limit, 
because mixing requires the decaying to zero of the cor- 
relations for all pairs of functions. Instead, a decay to 
positive values is observed if a suitable choice is made for 
the functions to be tested for autocorrelation. Actually, 
instead of considering the energies Ej of the single modes, 
we considered the energies of packets of n nearby modes, 
with n proportional to TV, precisely, the TVq quantities 

"(^■+1) TV + 1 

£j = ^2 Ek , where n = , TVq = 511 , 

I , 1 + 1 

for j = 0, 1, . . . , TVq — 1. In Fig. 4, right panel, the 
analog of c* for the quantities £j is plotted versus J/Nq, 
for the same numbers TV and K as in the left panel. It 
is true that the different curves do not superpose, and 
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FIG. 4: The "asymptotic values" c*{j) of the normahzed autocorrelation functions for the mode-energies Ej (left panel), and 
for the energies £j of packets of n nearby modes with n proportional to A'^ (right panel) , for an increasing number A'' of particles. 
The value of e is the same as in Fig. [2] The asymptotic value c*(j) is plotted in the left panel versus j/N for j = 1, . . . , N , and 
in the right panel versus j/No for j = 0, . . . , A^o — 1, with A^o = 511. In the inset of the right panel, for each of the four series 
of data the corresponding "moving averages" (over eleven points) are plotted. 



that a certain decreasing is observed, especially in pass- 
ing from TV = 511 to A*" = 1023. But for larger values 
of N the results seem to indicate that a nontrivial limit 
curve is being approached. This is better illustrated in 
the inset, where, in order to improve the readability of 
the graphs, the data were smoothed out by a standard 
"moving averaging" with eleven points. In our opinion, 
the results suggest that the FPU-type phenomenon dis- 
cussed here may persist in the thermodynamic limit, for 



a one-dimensional chain. And this, not for very special 
initial data, but in a global sense involving an averaging 
over all initial data, in a microcanonical setting. 

It may be of interest in this connection to recall that 
analytical perturbative estimates for a "freezing" in the 
thermodynamic limit, much in the same spirit of the 
present paper, i.e, by averaging initial data over an in- 
variant measure, were obtained very recently (see p7|). 
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